Driven flux-line lattices in the presence of weak random columnar disorder: 
Finite-temperature behavior and dynamical melting of moving Bose glass 
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We use 3D numerical simulations to explore the phase diagram of driven flux line lattices in 
presence of weak random columnar disorder at finite temperature and high driving force. We show 
that the moving Bose glass phase exists in a large range of temperature, up to its melting into a 
moving vortex liquid. It is also remarkably stable upon increasing velocity : the dynamical transition 
to the correlated moving glass expected at a critical velocity is not found at any velocity accessible to 
our simulations. Furthermore, we show the existence of an effective static tin roof pinning potential 
in the direction transverse to motion, which originates from both the transverse periodicity of the 
moving lattice and the localization effect due to correlated disorder. Using a simple model of a 
single elastic line in such a periodic potential, we obtain a good description of the transverse field 
penetration at surfaces as a function of thickness in the moving Bose glass phase. 
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Periodic structures driven on a random substrate, 
such as vortex lattices in type II superconductors, ex- 
hibit a rich variety of phases controlled by the inter- 
play between elasticity, disorder, temperature, and driv- 
ing forcei^^. At large enough velocity, where the ef- 
ficiency of the quenched disorder is reduced leading to 
dynamical reordering^, several moving glass phases were 
theoretically predicted^. In particular, considering elas- 
tic deformations for weak disorder or large velocity in 
d = 3, a topologically ordered phase is predicted^, 
namely the moving Bragg glass (MBG). The extension 
to correlated disorder led to the prediction of the mov- 
ing Bose glass (MBoG)£, characterized by the dynamical 
transverse Meissner effect (DTME), i.e. the tilt response 
to transverse field vanishes below a critical value. At 
finite temperature, renormalization group calculations^ 
show a transition at a critical velocity v c from the MBoG 
to a very high velocity glassy phase in which the DTME 
vanishes, namely the correlated moving glass (CMG). 
Clear evidence of the MBoG was found in numerical sim- 
ulations at T = 8 . However, a complete theory of elastic 
medium at high velocity is lacking, and the stability of 
moving elastic phases such as MBG or MBoG at the ther- 
modynamic limit is still debated. 

In this paper, we perform 3D molecular dynamics sim- 
ulations of superconductor vortices with weak random 
columnar pinning. We show in details the existence of 
the MBoG at finite temperature. It appears so stable 
that the expected dynamical transition to the CMG is 
not found. Furthermore, we find the existence of an ef- 
fective pinning potential that is z independent [z being 
the direction of the columnar pins) and periodic in the 
direction transverse to motion. Such effective pinning po- 
tential appears as an additional signature of the MBoG. 
Consequently, we extend to finite temperatures a model 
of a single elastic line into a tin roof potential 8 , which 
captures the DTME property of MBoG and yields quan- 
titative understanding of finite thickness effects. Finally, 
the dynamical melting of MBoG is studied as tempera- 
ture is increased. 



Following Ref4, we model a stack of N z Josephson- 
coupled parallel superconducting planes of thickness d 
with interlayer spacing s. Each layer in the (x,y) plane 
contains N v pancake vortices interacting with N p colum- 
nar pins parallel to the z direction. The overdamped 
equation of motion of a pancake i at position ri reads 

dr- 

= -YyiU vv {Pii,zn) - Yy^ vp ^v) 

+ F L +F tM (z) +Fi th (t) 

where P ij and z\j are the components of ry = rj — rj 
in cylindrical coordinates, pi P is the in-plane distance 
between the pancake i and a pinning site in the same 
layer at r p , and Vi is the 2D gradient operator act- 
ing in the (x, y) plane. r\ is the viscosity coefficient, 
F L = F L x is the Lorentz driving force due to an applied 
current, Fi is the thermal noise with (F* h a ) = and 
{Fl h a {t)Ff l {t')) = 2r 1 kBT8 ij 5 a p5{t-t') where a,(3 = x,y 
and fes = 1 is the Boltzmann constant. F tllt (z) is the 
surface force due to the field tilting away from the z axis 
in the y direction. This force acts as a torque on each 
flux line, i.e. F tm (z = 0) = -F tm {z = N z s) = F mt y 
and F ult (z) = for pancakes in the bulk. The tilt- 
ing force modulus is defined by F ult = e 2 cf)oHy/ '4ir = 

8ire 2 e \l b H y /y/3alH z , where e = (</>o/47rA a fc) 2 , ^ab is 
the in-plane magnetic penetration depth, ao is the av- 
erage vortex distance, H y is the transverse field compo- 
nent, and e is the anisotropy parameter. The intra-plane 
vortex- vortex repulsive interaction is given by a modified 
Bessel function U vv (pij) — 2e dK Q {pij/X a b)- The inter- 
plane attractive interaction between pancakes in adjacent 



layers of altitude z and [z + s) reads U vv (pi 



s) 



(2se /7r)[l + ln{\ ab /s))[{pij/2r g )' 2 - 1] for p y < 2r g and 
U^ipi^Zij = s) = (2se /7r)[l + ln(\ ab /s)][ Pij /r g - 2] 
otherwise ; in this expression r g = £ Q &/e, where £ a f, 
is the in-plane coherence length. This pairwise inter- 
action results from both electromagnetic and joseph- 
son coupling^. Finally, the attractive pinning poten- 
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tial is given by U vp (pi P ) = —aA p e (p*p/ r v> , where 
A p = (e d/2)ln[l + (-Rp/2^ h )}iS and a is a tunable pa- 
rameter. We consider periodic boundary conditions of 
(L x ,L y ) sizes in the (x,y) plane while open boundaries 
are taken in the z direction. Molecular dynamics simula- 
tion is used for N v vortex lines in a rectangular basic cell 
(L x ,L y ) = (5, 6V3/2)U o6 , with t = 1,2, and for a num- 
ber of layers varying from N z — 19 to N z = 1999. All 
details about our method for computing the Bessel po- 
tential with periodic conditions can be found in Rcf. 11 . 
The number of columnar pins is set to N p = N v . We 
consider the London limit k = \ a b/£,ab — 90, with an 
average vortex distance ao = X a b, and d = 2.83 10~ 3 A a fc, 
s = 8.33 lCT 3 A a6 , R p = 0.22 A a b, e = 0.01, r? = 1. 
We choose the tunable pinning parameter a — 1/25 so 
that the maximum pinning force is F^ ax ~ Fq/5 where 
Fq = 2eod/X a b is the unit force defined by the Bessel in- 
teraction. All the parameters values are identical to our 
T = previous study® so that a direct comparison is 
possible. 

The driving force F L is chosen high enough to obtain 
a fully elastic flow at T = 0. At low temperature, the 
vortex flow remains elastic with no dislocations. The 
rough static channels observed at T — persist, except 
that they are broadened by thermal fluctuations. Start- 
ing with such a lattice moving in the x direction and the 
magnetic field H along the z axis, we slowly vary the 
y component of H (H z being fixed) in order to obtain 
the transverse induction response. Fig. QJi shows the flux 
line inclination tan 0b = B y /B z averaged over time ver- 
sus the magnetic field inclination tan Oh = F[ y /H z for 
several thicknesses N z . The response is linear at low an- 
gle, with a finite slope which decreases when thickness is 
increased and eventually vanishes in the N z — > oo limit, 
as shown in Fig. [TJa. Such finite slope is explained in 
the insert of Fig. [IJl : the lines are curved at their ex- 
tremities by the tilt force while in the bulk they remain 
aligned with the z axis, i.e. the transverse field only 
penetrates the sample near its surfaces, resulting in the 
partial screening of H y . This supports that the finite re- 
sponse at low angle is a surface effect. Fig.QJi also shows 
that at a critical transverse field, the line inclination ex- 
periences a jump associated with an angular depinning 
transition of the vortex lines. Above this transition, the 
lines display a kink structure in the yz planes (see Fig. 
3 of Ref4). dz/dy is almost independent of the line for 
a given y and periodic in the y direction, indicating that 
all the lines experience the same effective pinning land- 
scape, which has the invariance of the columnar pins and 
the transverse periodicity of the channels. These results, 
very similar to those obtained at T = 0, strongly suggest 
the existence of the MBoG phase at finite temperature. 
Since a finite transverse magnetic response at low angle 
is expected in the CMG£, a careful study of finite size 
effects in the z direction is crucial to discrimate MBoG 
from CMC In order to understand quantitatively the 
influence of thickness, we study a model introduced in 
Ref. 12 and used in Rcf. 8 to describe the tilt of the lines 
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Figure 1: (color online) (a) Average vortex line inclination 
versus field inclination at T = 10~ 4 for v = 10~ 2 , N v = 
30 and several thicknesses (circles) compared with disorder 
free linear response (squares). Insert : average shape of a 
line at low inclination for N z = 149 and T = 10~ 10 (thick 
brown line), and sinh fit (thin dashed blue line), (b) Slope 
d(tan#.e)/d(tan#jf) at low angle versus N z for T = 10 -10 
(red circles) and T = 10 -4 (green triangles). The solid blue 

line is the tanh (^j^j /L fit, the dashed green line is the iV" 1 
dependence at large N z . 



at T = 0, that we phenomenologically apply to quanti- 
ties averaged over time in the finite temperature case. It 
is a mean field approach describing the angular response 
of the vortex lattice in terms of a single elastic line put 
in the effective pinning potential V{y) discussed above. 
The energy E of a line of length L = N z s is given by 



E{u) 



Aim 



V{u) + f (u(L) - u(0)) 



where u(z) is the one-dimensional displacement field in 
the y direction, c = e 2 e is the elastic constant and 
/ oc H y /H z is a surface force. This expression of 
the energy doesn't take explicitly into account ther- 
mal fluctuations, however line elasticity and effective 
potential can depend on temperature. Minimizing E 
with respect to u(z) while expanding V quadratically 
near a minimum leads to the following solution u{z) 

for a line : u{z) = ^- sinh ( z z J 2 ^j / cosh C^^j, where 

zq = \J c (^^") u _ characterizes the thickness of the re- 
gion where the transverse field penetrates the sample 
near the surface. The average line inclination is given 

by tan6» B = \u(L) — w(0)] /L = ^tanh^). The 

above equation for u(z) is observed to accurately fit the 
data, as shown in the insert of Fig. QJi. z calculated 
from this fit is found to be independent of both N z and 
6h, which is consistent with the observation that at low 
Oh (i.e. when u(z) is small enough for the quadratic 
expansion to be accurate), tan#s is a linear function of 

tanOn (see Fig. QJi) and tan#B oc tanh j /L (see 
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Fig. |T|b) which is verified up to large thicknesses com- 
pared with the penetration length. In the large thickness 

limit L zq, tanh (^j^j ~ 1 and the inclination scales 

as L^ 1 , i.e. as iV" 1 . From this we conclude to a true 
vanishing response to transverse field below a threshold 
in the infinite thickness limit, confirming that we are see- 
ing MBoG and not CMC Finally, our results support the 
existence of the effective pinning potential at T > and 
reinforce the interest of the simple mean field model — 
extended to finite temperatures provided that simulated 
quantities are averaged over thermal fluctuations — as 
a tool to describe DTME, including finite size effects. 
These behaviors have been observed in a wide range of 
velocities, from the appearance of the MBoG to the high- 
est velocities we can simulate because of duration issues 
(over more than 4 orders of magnitude at T — 10 -7 , over 
3 orders of magnitude at T — 10 -4 ), i.e. we don't see 
the expected dynamical transition to CMG 7 -. Since the 
critical velocity v c at which the transition to CMG is ex- 
pected behaves like 1/L C oc N p a 2 where L c is the static 
Larkin length, a similar analysis has been performed for 
weaker pinning strengths. No evidence of CMG phase 
has been found in that case either. 
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Figure 2: (color online) Cv = ¥ yV with Hb = 1 for 
N z = 49 (blue squares) and slope of the transverse induction 
response d(tan 6b)/ d(tan 6h ) at low angle versus temperature 
for N g — 19, 29, 49 (red circles) compared with the pinning 
free case (green triangles), both for Nv = 30 and v = 10 -2 . 

We now fix the velocity (v = 10 -2 ) and study the ef- 
fect of temperature. At zero tilt {Oh = 0°), we compute 
C v = {{El) ~ {E t ) 2 )/k B T 2 (Ref. 13) where E t is the 
total interaction energy {Cy would be the specific heat 
at thermodynamic equilibrium). Together with Cy, we 
plot in Fig. [5] the slope d(tan b) / 'd(tan 6 h) of the lin- 
ear region of Fig. |TJl at low angle. At low temperature 
d(tan#s)/<i(tan#ff) is reduced compared with the pin- 
ning free case and decreases when thickness is increased, 
illustrating the DTME property as seen previously. It 
experiences a jump around T = 2.5 x 1CT 4 while the 
thickness dependence vanishes, what we interpret as the 
loss of DTME. Concomitantly, Cy exhibits a sharp peak 
suggesting a dynamical phase transition. To elucidate 
the nature of the phase obtained once MBoG has disa- 
peared, we compute in the xy plane the pair correlation 




Figure 3: (color online) (a) Cv (blue squares), height of 
the correlation function first neighbor peaks (red circles) 
and background level (green triangles) versus temperature 
for N v = 120, N z = 19, v = 10" 2 . (b) Amplitude of 
the in plane correlation function along the x axis just be- 
low (T = 2.1 x 10" 4 , red line) and above (T = 2.2 x 10" 4 , 
dashed blue line) the transition. 

function g{r) = (p(r' + r)p(r')) (Ref. 13), where p{r) 
is the vortex density, and the mean square displacement 
B x {t) = ([r(t) - r(0)] 2 ) (Ref. 13), both at zero tilt. At 
low temperature, the correlation function displays a tri- 
angular lattice of peaks, in agreement with the expec- 
tation of a quasi ordered lattice for an elastic moving 
glass, while above a critical temperature we find circu- 
lar rings and short range order, signature of a disordered 
isotropic phase. In order to compare the critical temper- 
ature obtained from the correlation function and the one 
obtained from Cy , we calculate the height of the peaks 
corresponding to the six first neighbors and the back- 
ground level. As shown in Fig. [3^, the collapse of these 
two quantities — which indicates the loss of the sixfold 
symmetry — occurs at the same temperature Tc as the 
peak in Cy. In Fig. [3)d we plot the correlation function 
along the x axis just below and above Tc, pointing out 
the change in the order range. This picture is confirmed 
by the study of the mean square displacement. Below 
the critical temperature B\{t) is bounded, while above 
it grows linearly indicating a diffusive wandering as ex- 
pected in a liquid. However, the diffusion coefficient is 
found to be a nonlinear function of the temperature, sug- 
gesting that the motion is more complex than a classic 
random walk. To sum up, at a critical temperature Tc 
we simultaneously see the vanishing of both the DTME 
and the lattice order, while the vortex displacements go 
from bounded to diffusive. We conclude that the transi- 
tion observed is the dynamical melting of the MBoG into 
a moving vortex liquid (MVL). Finally, we take a look at 
the correlations along the vortex lines in order to high- 
light the influence of the effective pinning potential V . 
We plot in Fig. 0^ the mean square displacement in the 
transverse direction within the lines Bi{zij) = {{yi—yj) 2 ) 
(Ref. 13,14) where i and j are two pancakes belonging to 
the same flux line and z%j is the distance between the two 
layers they belong to. At low temperature B2 is bounded, 
each line being pinned in a minimum of V. This pinning 
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Figure 4: (color online) (a) B2{z%j) = — j/j ) 2 ) (Zij is ex- 
pressed in units of layer spacing s) for T from 10 -4 to 1CP 3 , 
N v = 30, V z = 149, v = 10" 2 . Insert : B 2 (zy) versus T for 
Zij = 10s (green circles) and Zij = 100s (red squares), (b) 
From top to bottom : total energy, amplitude of the correla- 
tion function at the first neighbor peaks and i?2(L/2) versus 
time close to melting for N v = 120, N z = 29, v = 10~ 2 . 



effect weakens as temperature is increased, and eventu- 
ally vanishes at T = Tq- The insert of Fig. [4Jl clearly 
displays the jump in B2(T) associated with this transi- 
tion. The high temperature behavior can be understood 
by assuming that only thermal fluctuations and line co- 
hesion are relevant : the length of each bond between 
two neighboring pancakes is then an independent ran- 
dom variable, and the line configuration is analog to a 
random walk in which z would be the time, leading to 
the observed linear growth. We conclude that the loss 
of localization along the z axis is a manifestation of the 
disappearance of the effective pinning potential. To be 
sure that this transition and the in-plane melting are two 



aspects of the same phenomenon, we monitor the time 
evolution of the quantities indicating the transition. Be- 
cause of finite size, close enough to the transition the 
system hesitates between different phases, and continu- 
ously switches from one to another. We can check in 
Fig. 0Jd that the three indicators of the transition plot- 
ted versus time (total energy, amplitude of correlation 
function at first neighbors peaks and B2{L/2)) oscillate 
between two states, which we identify respectively with 
MBoG (lower energy) and MVL (highest energy), and 
that in-plane and out-of-plane quantities jump simulta- 
neously when the system goes back and forth from one to 
the other. The effective pinning potential thus persists 
in MBoG whatever the temperature, and only vanishes 
when melting occurs. 

To conclude, we find evidence of the MBoG phase 
at finite temperature, exhibiting DTME below a critical 
transverse field. It is stable in a wide range of tempera- 
ture and velocity, respectively up to the melting tempera- 
ture and to the highest velocities we are able to simulate. 
Weaker pinning, which should lower the critical velocity 
v c , has also been studied, but no evidence of the CMG 
phase has been found in that case either. A reduced den- 
sity of pinning centers should also lower v c , and could be a 
direction to further investigate the existence of CMG. We 
also predict the existence of an effective transverse static 
tin roof pinning potential in the MBoG phase. Since the 
CMG doesn't exhibit DTME, we expect this potential to 
vanish in the CMG as it does in the MVL. 
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